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PREFACE 

The Coasul Zone Color Scanner (CZCS) on board the 
Nimbus 7 satellite was designed to measure surface radiance upwelled 
from the ocean in 6 spectral bands. The CZCS spectrometer obtains 
its information from a rotating mirror and is timed to collect data 
when the mirror views the earth surface between ca. 40° to the left 
and right of the subsatellite track. Each scan is divided into 1968 pic- 
ture elements, pixels, of 0.04 degrees scan each. In order to avoid 
direct reflected sun glint, the rotating mirror shaft can be tilted so that 
scans cross the subsatellite track up to 20° forward or aft of the point 
directly beneath the satellite. The CZCS is the first satellite borne 
instrument to have this tilted scan capability and therefore poses some 
new problems in locating the earth surface position of viewed pixels. 
This report resutes results of analyses performed by several members 
of the Nimbus Experiment Team for the CZCS, using a consistent 
coordinate system and terminology. 

This work was supported by the National Aeronautics and Space 
Administration, Contract No. N \S 5-26249, and the US E)epartment 
of Commerce, National Oceanic and Atmospheric Administration, 
Grant No. NA80AA-D-00007. 
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THE CZCS GEOLOCATION ALGORITHMS 


The location on e^h of a pixel viewed at a ceruin tilt and scan angle is a complex problem in CZCS view- 
ing geometry and spherical trigonometry. Presented here is the derivation of the Visibility Laboratory algorithm 
for CZCS pixel location. It w originally developed by Wayne Wilson and Raymond Smith of the Visibility 
Laboratory , in collaboration with Jim Mueller of the Naval Post Graduate School, Monterey and Howard Gordon 
of the University of Miami. 

SATELLITE GEOMETRY 

Rg. 1 defines the coordinate system used relative to the satellite. Note: angles are symbolized by greek 
letters, vectors are .symbolized by lower case letters superscripted with a bar (a), and rotation operators, (matri- 
cies) are symbolized with capital letters. 



Rg. 1 


DEFINITIONS 

X axis The x axis is along the direction of flight of the satellite, positive towards the radiometer sensor. 

z axis The z axis connects the center of the scan mirror and the center of the earth, positive towards the earth 

and orthogonal to x axis at the point of intersection. 

y axis The y axis is orthogonal to the x and z axes. 

a a is the vector through the shaft of the scan mirror (the mirror rotates about this vector). Vector a 

rotates about the y axis in the x-z plane defining the tilt angle r. 

n n is a unit vector orthogonal to the surface of the scan mirror. 

0 0 is a unit vector along the x axis (oriented toward the sensor). 

V V is j unit vector from the direction of the pixel to the scan mirror. 

Note that vectors 6, n, and v lie in the same plane. 

9m 9„ is the angle between o and n. 

r T is the tilt angle of the scan mirror shaft. (One half the "tilt" given in the subcommutated data ) 

n Tj (not shown) is the mirror shaft rotation angle. 
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The tilt and scan angle are known parameters. From these we would like to calculate the nadir angle 0 from 
the z axis to vector v, and the azimuth angle ^ from the x axis to the projection of v on the x*y plane (Fig. 2). 



Fig. 2 


For convenience in the derivations to follow it is desirable to rotate the coordinate system above to a 
^natural* coordinate system such that z„ - x, x„ - -z. and y„ - y (Fig. 3). In this system the angles describing 
V, (#„ and <t>„) can be more easily calculated from t and 17 After v is described in terms of and < 6 n the coordi- 
nates will be rotated to the "satellite" system and 6 and <t> derived. 



Recall that vector h is the normal to the scan mirror and therefore is dependant only on tilt and scan for its 
orienution. The initial position of vector n at zero tilt and zero scan can be defined as lying in the x„-z„ plane at 
an angle 45* clockwise rotation from the z„ ax>s, (Fig. 4). 



Fig. 4 
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Since n is a unit vector, its x„, Zn coordinates may be calculated: 


no- 


X„ 


COStf„ 

-1 / V2 


-a 

Xn 


0 



0 

Zn 


sin«„ 

1 1 / V2 


a 


The CZCS mirror tilts about the y„ axis, (assumed to be through the surface of the mirror). The tilt angle r 
is defined positive when n rotates in a counter clockwise direction. An Eulerian rotational operator T can be 
defined such that it acts on n to cause a counter clockwise rotation about y^. 


T 


COST 0 sin T 
0 1 0 
-sin T 0 COS T 


The mirror scans about the mirror shaft; at zero till this shaft lies along the Zn axis. An Eulerian rotational 
operator S can be described to rotate n around the z„ axis. When facing the negative z„ axis direction, this rota- 
tion is positive in a clockwise direction. 


S 


COS 17 sin 7) 0! 

—sin 7) cos 7) 0 
I 0 0 li 


We would like to combine the operators into a single operator. Normally one considers the CZCS mirror to 
scan about a tilted mirror shaft, however if we operate T on n the mirror shaft (fixed in relation to n) also rotates 
out of alignment with Zn and the use of S to scan about z^ becomes invalid. It is mathematically valid to scan n 
about z„ first then tilt it about y„ second. Vector n will have the same final coordinates and the simpler Eulerian 
operators can be used. 


TS - T (S) 


cos r cos 7) 
—sin 7 f 

—sin T cos 7j 


cos T sin 7 f 

cos 7) 

—sin T sin rj 


sin T 
0 

cos r, 


Therefore for any t and tf. 


n' ^ 

Xn 

Vn 

-TS 

-a 

0 


a (sin T - cos T 
a sin 17 

cos T)) 


Zn 


a 


a (cos T + sin T 

cos T)) 


We now have the coordinates ^or the mirror normal vector n after scan and tilt and can proceea to derive 
the coordinates of vector v from the direction of the pixel. Recall that vector o is a unit vector lying along the z„ 
axis in the direaion of the CZCS sensor. 


0 - 


0 

0 

ll 


Vector 0 is separated from vector n by the angle 
the angle of reflection, vector v is also separated from n by 


Because the angle of incidence to a mirror is equal to 
angle 9 ^ and lies in the ssme plane as vectors 0 and n. 
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Since vector v is directed into the mirror and vector o is directed away from the mirror, if vector v is subtracted 
from vector o the resulting vector will lie along n with a length of 2 cos 9m . (Fig. 5). 

o - V - n 2 cos 9m 
V - 0 - n 2 cos 9m 


V 


0 - 2 cos dm • a (sin t ~ cos t cos t)) 

0 — 2 cos dm a sin n 

1 - 2 cos d„, • a (cos t sin r cos 17 ) 


With the coordinates of v, angle can now be calculated. 


Xn sin n 

x„ sin r - cos T cos Tj 


Note that computational machines return arctangent functions in the range of ±90®. Therefore for normal scan 
angles 180* must be added to the returned value of (Fig. 6 ). Therefore: 


“ 180® ± arctan(sin tj / (sinr - cos t cos tj) ) 


Since lo| - |nl - 1, 


n • o ■■ |6l |nl cos dm - cos dm • a(cos t ± sin r cos n) . 
d„ - 2dm 
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- 2 arccos (a(cos t + sin r cos tj) ) 



Vector V is now described in the "normal" coordinate system in terms of tb^ angles and To rotate 
back to the "satellite" coordinate system recall tiiat; 


"satellite" "normal" 


X 


sin 9 QOS 


Zn 


COS 9n 

y 

- 

sin 9 sin <f> 

* 

Yn 

- 

sin 9n sin 

z 


cos 9 




-sin 9n cos <6n 


Therefore: 


9 • arccos (-sin cos 


and 


- arctan(y/x) - arctan(tan 9^ sin <tn) • 


Due to the quadrant ambiguity of arctangent function, for negative tilts, 180® may need to be added to to 
correct 

PITCH - ROLL - YAW 

The "satellite" coordinate system as defined above is fixed relative to the space craft. However, the satellite 
may be rotated due to pitch, roll, and/or yaw relative to an "earth" coordinate system. The origin of the "earth" 
system is the same as for the "satellite" system, the center of the satellite. The x direction is the instantaneous 
tangent to the flight path, the z direaion is towards the center of the earth, and y is orthogonal to x and z. The 
"satellite" system will oscilate about the "earth" system by small angles generally less than a degree. Fig. 7 f ows 
the directions of the pilch, roll, and yaw rotations on the "earth" system, these rotations define the orientation of 
the "satellite" system in Fig. 2. 
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Fig. 7. j 


At this point we have vector v defined by angles <t> and 9 relative to the "satellite" coordinate system. The x, 
y, and z coordinates of v may be calculated as in Eq.d) above: 

X - sin ^ cos 0 
y — sin B sin <t> 
z - cos d . 

We now need to create a set of operators that will rotate the "satellite" x, y, z coordinates to a set of "earth" 
x\ y\ t! coordinates for v. The order of correction is significant. The xyz coordinates must first be pitched about 
the y axis, then rolled about the new x axis and finally yawed about the again new z axis. First the pitch operator: 

I cos p 0 sin p 

0 10 , 

-sin p 0 cos p 

where p equals the pitch angle. Next the roll operator: 


1 

0 

0 

0 

cos r 

-sin r 

0 

sin r 

cos r 


Finally the yaw operator: 


cos y 

-sin y 

0) 

sin y 

cos y 

0 

0 

0 

1 


These may be combined by first multiplying R on P and then multiplying Y on RP yielding YRP 


cos y cos p + sin y sin r sin p -sin y cos p + cos y sin r sin p cos r sin p 
YRP • sin y cos r cos y cos r -sin r 

-cos y sin p 4* sin y sin r cos p -sin y sin p + cos y sin r cos p cos r cos p 
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Sine? we are dealing with small angles, considerable compuution time can be saved by assuming the cosine 
of f; p, or y — 1, (cos 1® — 0.9998), If working in the radians mode one can assume that the sine of r, p, or y 
equals r, p, or y. Therefore, YRP becomes: 


YRP 


1 + yrp ~y + rp p 
p 1 -r 

-p + yr -yp + r 1 


If working in the degrees mode, first precompute sin r, sin p, and sin y. 

Now we may multiply YRP on a set of "satellite" coordinates and calculate a set of "earth" coordinates for v. 
For clarity let the satellite x, y, z coordinates be expressed as u, v, w to distinguish them from the angles r, p, and 
y. Then: 


x' - (I -f yrp)u + (-y + rp)v + pw 
y' — yu + V + rw 
z' “ (“P + yr)u + (-yp + r)v + w . 


Notice that if roll, pitch, and yaw arc zero, x' ■■ u, y' v, and z' — w as we would expect. 
Finally, we can calculate 9' and for the "earth" system: 


9* - arccos (z') 


- arctan (y'/x') . 

Due to the quadrant ambiguity of the arctangent function, for negative tilts, 180® may need to be added to <^' to 
correct 

EARTH LOCATION 

If the latitude and longitude of a geographical position, and the distance and true bearing to a second geo- 
graphical position are known, the latitude and longitude of the second geographical position can be calculated. 
Distance and true bearing from the satellite ground point to the pixel can be derived from 9 and or 9' and if 
correcting for roll, pitch, and yaw. 

Given the orbital altitude of the satellite, h, for example 952km, and the nadir angle 9, the surface distance, 
6, in degrees great circle arc is: 


8 - arcsin(R sin 9) - 9 


where 


R- (r + h)/r- 1.1494, 


and 


r - earth mean radius - 6371.2km . 


Tlie inclination of the CZCS orbital trackline is a function of latitude, and is 9.28® (west of north) at the 
equator. For any latitude, the inclination, a, is: 


a - arcsin(sin 9.28/cos lat) . 


True betting, /3, then is ^ - a, if /3 is negative, ~ 360 + j8. 

Pn 



Fig. 8 represents a general spherical triangle. 


Pmnts: 

G.P. G.P. is the ground point, subsatellite point, of the satellite. 

Pn Pn is the north pole. 

Pi Pi is the pixel. 

Angles: 

/9 j9 is the true bearing from the satellite G.P. to the pixel. 

^on Alon is the difference in longitudes between points G.P. and F. 

S 8 is the great circle distance between G.P. and F. 

Colato Colato is the co-latitude of the G.P. (Colat ” 90* • lat) 

Colati G)lat| is the co-latitude of the pixel. 

If the latitude and longitude of the G.P. are known then the latitude and longitude of the pixel may be cal- 
culated directly. From the cosine rule: 

lat| - arcsin(sin lato cos 8 -I- cos lato sin 8 cos /3) . (2) 

A lat or Ion subscripted with a 0 (lato) refers to the G.P. and those subscripted with a I (lat|) refer to the pixel. 
Notice that the sine of the colatitude equals the cosine of the latitude, etc. From the sine rule 

dion - arcsin(sin sin 8/cos lat|) . (3) 


Then, 


looi “ lono - Alon 


These calculations work for all latitudes north (-!■) or south (-) to ±80.72* and all longitudes, (4 west, - east), 
for the north bound satellite track. 


8 


If the latitude and longitude of the pixel are known, then the calculation of the latitude of the G.P. is a 
complex quadratic equation since the true bearing, /3, is a function of the unknown lato. From above: 


Let 


^ « d - arcsinfein 9.28/cos lat<>) 


Then 


I « sin 9.28/cos lato • 


(4) 


/3 - d - arcsin I, 


and we may substitute this in Eq. (2) above. 

lati arcsin (sin lato cos 5 + cos lato sin 8 cos(d - arcsin I) ) . 


The term containing I may be rewritten: 


cos(d ~ arcsin I) — cos d cos arcsin I + sin d sin arcsin I 


Note that 


(5) 


cos arcsin 


i-vr^ 


Therefore Eq. (S) equals: 


•• cos d Vl — P + I sin d 


“ cos d VI - sin* 9.28 / cos^ lato + sin d sin 9.28/cos lato 


- lato - sinJ 9.28 + iiUJUiniiS 

cos lato ^ cos lato 


Substituting this last back into Eq. (2) the lato terms cancel yielding: 

sin lat| - sin lato cos 8 + sin d (cos d >/cos^ l«to - sin^ 9.28 + sin d sin 9.28) . 

This can be rearranged yielding: 

sin lato cos 8 + sin 8 cos d VI - sin* lato - sin' 9.28 - sin lat| - sin 8 sin d sin 9.28 . (6) 
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Ui 


X - sin lato 

a - sin lati - sin 5 sin 0 sin 9.28 
b - cos 5 
c « sin d cos ^ 

d - 1 - sin^ 9.28 « cos^ 9.2^ . 
Equation (6) becomes: 


bx + c Vd - x^ - a . 


Rearranging and removing the radical yields: 

x^b^ -h c^) — x(2ab) + (a^ - cM) • 0 


Thus: 


ab ± Va^ - (b^ c^) (a^ - c^ d) 

b^ 


The value under the radical is chosen to be positive if tilt, (r), is negative and negative if tilt is positive. Then 

lato * arcsin x . 

Note: It is simplest to precompute the values of a, b, c, and d before attempting to calculate x. 

Beta is now calculated as in Eq. (4), and Alon is calculated as in Eq.O). Lono then is: 

lono • loni Alon . 


TABLE OF USEFUL INFORMATION 

Mean satellite orbital period (Nov- Dec 1980) - 104.07 minutes 

Mean orbiul velocity (great circle arc) — 5.765 x 10'^ degrees/ second 

Mean satellite orbital altitude (Nov-Dec 1980) * 952 kilometers 

Satellite scan rate * 123.75 milliseconds/ scan line 

Trackline inclination angle « 9.28 degrees west of north 

Mean earth radius 6371.2 kilometera 

Earth rotational velocity - 4.1667 x 10"^ dcgrces/second 
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PROGBAM LISTINGS 

The foUowing are listings of two sets of programs using the CZCS geolocation al^rithms. First is a HPL 
program for the H.P. 9825 calculator. This program is primarily a set of subroutines which may be used for spe- 
cial stand akme applications. The calling routines, accessed by special function keys, provide basic calculations and 
may be easily modified by for-next loops to provide for example, trackline latitudes and longitudes. 

Next, is a set of Fortran programs used in the Visibility Laboratory’s CZCS atmospheric correction algo- 
rithm. These subroutines, "CZSUBr - *CZSUB5* and *C2 !SI}ba'', provide all the basic algorithm calculations. 
Also included is "CZNAV*, a utility program similar to the HPL program, which calls the various CZ Subroutines. 
All the Fortran programs are written to compile on a Prime SSO Computer. 'C25NAV* a written in Foruran 77. 
The CZ Subroutines have been compiled using Fortran 77 but %vere written using Fortran 66 protocol. 

Or "CZNAV: 

1: "Stand alone CZCS tracking program": 

2: "09 January 1981, visLab, Nolten": 

3: aim L(6] ,H(2] ,C$[32] 

4: getK "CZKEY* 

5: 6*F;ent "Enter output device. ",P;gsb "Notes" 

6; gto "START" 

7: 

8; "- — -" :f or K»1 to 35;wtb F," ";next K;wtb F," ",13,10;ret 

9: 

10: "dm" : 

11: "returns angles in degrees. minutes oecimal minutes”: 

12: prnd(pl,-6)«pl;ret int (pl)-*-.6trc( pi) 

13: "dd": 

14: "returns angles in decimal degrees": 

15: ret int (pi) -fire (pi)/. 6 
16: 

17: "dist": 

18: "pi * nadir angle; returns p2 » distance": 

19: asn( ( (952+6371)/6371) sin (pi) ) *p2 ;p2-pl*p2 

20: ret p2 

21: 

22: "pix": 

23: ”pl«tilc, p2*pixel i, returns: p3Bazimutn, p4«nadir ang, p5*snatt ang”: 
24: -39.04-*-. U4p2«p5 

25: 2acs( (cos (pi) -fsin(pl) cos ( p5) ) /V2) ->p6 

26 : 18U-*-atn(sin ( p5)/ ( sin(pl)-cos( pl)cos(pS) ) ) -*p7 

27 : " tneta" :acs(-sin (p6 )cos(p7 ) ) -»p4 

28: if p6*90 or p6*270;sgn (p5) *2*l8U-9Usgn( p5) ->p3; jmp 2 
29: "pni" :atn (sin(p6) sin(p7)/cos(p6) )-*p3;i£ pKO; 180-*-p3-*p3 
30: if not (m or X or Y);ret 
31: "Roll-Pitcn-Yaw" : 

32: "x,y, 2 “ :sin (p4) cos( p3) -*p8;sin(p4) sin(p3) -*p9 ;cos(p4) -*pl0 
33: "sin (R. P. Y ) " :sin( w)-»pll ; sin (x) -»pl2 ; sin ( y) -»pl3 
34: "x* " : (l-*-pllpi2pl3) p«*-(-pl3-t-pllpl2)p9-*pl2pi0-*pl4 
35: "y ’ " :pi3pb-*-p9-pilpi0-*pl5 

36 : "2 ' " : (-pl2-*-pllpl3) p8-* (-pl2pl3-*-pll) p9-*pl0-»pib 
37: "tneta* " :acs( pie) -*p4 

38: if pl4;atn (plS/pl4) -»p20 : if T<0; 18U-*p20->p20 
39: if not pl4: sgn (pl5) 9U-»p20 ; if KO ; 180-i-p20 ->p20 
40 : "pni ' " :p20*p3 
41: ret 
42: 

43: "pixdir": 

44: "entry: pi-tracx line azimuen, p2»G.P. lat: return: p3«true bearing": 

45: pl-aan(co8( 8u.72)/cos(p2) )*p3;if p3<0; 3b0-*p3-*p3 

46: ret p3 

47: 

48: "pixlat”: 

49: "entry: pl"G.P. lat, p2"Sfc dist, p3>true bearing; return: p4*pix lat": 


II 


a 


50 : 9 0-pl ♦pS ; acs ( cos (p5)cos(p2)+8in(p5)sin(p2)cos(p3))^p4;9 0-p4 *p4 

51: ret p4 

52: 

53: "pixlon": 

54: "entry: pl«Pix. lat, p2«G.P. Ion, p3«stc dist, p4«true bearing": 
55: "return: p5*pixel longitude": 

56 : p2-asn ( sin ( p4 ) s in ( p3 ) /s in ( 9 0-pl ) ) *p5 

57: ret p5 

56: 

59; "equa" : 

60: "entry:pl*GP lat, p2»GP Ion; return: p3»sec since, p4«deg since": 
61: "p5»scan line since, p6-longitude of equator": 

62: asn(sin(pl)/sin( 80. 72) ) *p4 ; p4/5. 765e-2*p3 

63: p 3/.12 375*p5;p2-asn(tan(pl)/tan(80. 72) ) -4. 1667e-3p3*p6;ret p6 
64: 

65: "scanli": 

66: "entry: pl*seconds, p2»aegrees; return: p3*scan lines": 

67; if p2t0;p2/S.765e-2*pl 
68: pl/.12375*p3;ret p3 
69: 

70: " second" : 

71; "entry; pl»scan lines, return: p2“seconds, p3»degrees"; 

72: .12375pl*p2;p2"5.765e-2*p3;ret p2 
73: 

74: "latlon": 

75: "entry; pl*new sec since equa, p2*lon of equa"; 

76: "return; p3«new lat, p4*new Ion": 

77; 5. 765e-2pl*p5;asn( sin (80. 72) sin (p5) ) •►p6 

78 : p2+asn(tan( p6)/tan( 80.72,' ) +4. 1667e-3pl*p4 ; p6-*p3; ret 

79; 

80: "GP": 

81; "pl-Aziinuth , p2>Oistance, p3»Lat-l; ret p4*Lat-0": 

82: sin (p3)-sin(p2) sin (pi) sin (9. 28) «rl 

83: cos(p2)*r2 

84: sin(p2) cos(pl) -»r3 

85; cos(9. 28) »2*r4 

86: rlr2*r5; r2»2+r3*2*r6 

87: r5*2-r6(rl»2-r3*2r4) *r7 

88: if r7>0;Vr7*r7 

89: if not (sgn (T) -*r8) ; l*rb 

90; asn( (r5-sgn(r8) r7) A6) *p4;ret 

91: 

92: "I ter” :Cxo 0 
93: 0-*I*0;1969*V 

94: "LI" : (UeV)/2-*P;cll ’pix' ( 'l/2,P, A,n,R) 

95: cll ‘GP* (A, 'ai8t'(N,D) ,L(31 ,L(1) ) 

96: cll 'pixair*(A,L(l),d) ; L ( 4 ) +asn(s in (B )sin (D) /cos (L (3 ) ) )*L(2) 

97: cll •equa’(Lll),L(21,S,0,M,K) 

9b: I*l*l;asp "Iter;",I;if I>20;999*U 
99: it abs(K-L*r0) <.00028 or R«999;gto ■♦■4 
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100: if r0<0;P«U 
101: if r0>0;P«V 
102: gto "LI" 

103: ret 
104: 

105: "START" : 

106: fnt c35," • ",fl5.4 

107: beep;dsp "Key Special Function to Continue" ;stp 
108: 

109: "5": 

110: "DegHin-DecOeg and ScanLine-Time conversions": 

111: pet "Conversion Codes" ;fmt clO," »",c4 

112: wrt 16 , "seconds", "S" , "scan 1 ines", "L" , "dec deg" , "DO" , "deg min", "DM" 

113: wrt 16,"U0IT","u";spc 2 

114: ent "Enter value and Identity code",C$ 

115: if. pos(cap(C$) ,"u");cll • ';gto "START" 

116: vaKCSX 

117: fmt l,10x,f7.2," SECONDS «",f6.0," SCAN LINES" 

118: fmt 2,10x,f7.0," SCAN LINES «",f7.2," SECONDS" 

119: fmt 3,10x,fl(i.4," DEGREES MINUTES >",fl0.4," DECIMAL DEGREES" 

120: fmt 4,10x,fl0.4," DECIMAL DEGREES -",flC.4," DEGREES MINUTES";fxd 1 
121: it pos(cap(C$) , "S" ) ; wrt F+. 1,C , 'scanli* (C) 

122: if pos(cap(C$) ,"L") jwrt F+. 2, C, ’second' (C) 

123: if pos(cap(C$) ,"M") ;wrt F+. 3,C , ’dd’ (C) 

124: if pos(cap(CS) ,"DO") ;wrt F+. 4 ,C , 'dm* ( C) 

125: gto -11 
126 : 

127: "6": 

128: "Enter Equa crossing Lon & Time, Sta. Lat-Lon, get sta pix vals 6 time": 
129: ent "Enter Equatorial crossing Long.”,L 
13U: ent "Time ct crossing (Hri.MMm)",Z 

131: ent "Enter Station Latitiude",L(3) , "Station Longitude", L(4J , "Tilt" ,T 
132: gsb "Iter" 

133: wtb F,10,10 

134: wrt F, "STATION LATITUDE" , 'dm ’( L ( 3J ), "STATION LONGITUDE" , 'dm' ( L (4 j ) 

135: wrt F, "EOUATORIAL LONGITUDE" , ’am’( L) , "TIME AT EQUATOR", 2 
136: wrt F, "TILT",T 

137: Cll ' ';if R«999;wrt F, "STATION NOT SEEN THIS PASS, R",999;gtO ♦4 

138: wrt F,"TIME OF STATION PASSAGE" ,’ dm ’(' dd '( 2 )♦' dd '( S/6000) ) 

139: wrt F, "STATION PIXEL NUMBER" , prnd (P ,0) 

140: wrt F, "SECONDS SINCE EQUATOR" , S ,"SC AN LINES SINCE EQUATOR" ,M 
141: gto "START" 

142: 

143: "7" : 

144: "Enter Pixel Lat-Lon, Pixel I, get G.P. Lat-Lon, & Equa aata": 

145: ent "Enter Pixel lat . " , L( 3) , "P ixel Ion. " , L( 4 ) , "P ixel » ", P, "Tilt" ,T 
14 e : cli 'pix' (T/?,P,A,N,R) ;cll 'aist'(N,D) 

147: cli 'Ct-' (A,D,M3| ,L(1)) 

148: cli 'pixoir'(A,L(ll,J);L(4|4’asn(sin(B)sin (D )/cos (L( 3| ) )*L( 2) 
i49: cii 'equa' (L[i| ,L(2) ,S,0,M,L) 
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150: 

151: 

152: 

153: 

154: 

155: 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 
18U 
181 
182 

183 

184 

185 

186 
18 7 
188 
188 

190 

191 

192 

193 

194 

195 

196 

197 

198 
198 


wtb F,10,10 

wrt P, “PIXEL LATITUDE", 'da* (L( 3] ), “PIXEL LONGITUDE" , *da* ( L[4] ) 
wrt P, “PIXEL #“,P,“TILT“,T 

Cll • ’>wct P, “G.P. LATITUDE", 'dm* (L(1J) 

wrt P,“G.P. LONGITUDE", 'da* (L (2]), “EQUATORIAL LONGITUDE" , 'da* (L) 
wrt P, "SECONDS SINCE EQUATOR" ,S ," SCAN LINES SINCE EQUATOR" ,N 
gto "START" 


" 8 " : 

"Enter Pixel Lat>Lon, X~Y corrds, second pix x-Y coords,*: 

"Get second pixel Lat-Lon" : 

ent "Enter 1st Pixel lat.*,L(3) ,”lst Pixel lon.*,L(4] ,*lst Pixel t*,P 

ent "1st Pixel Scan line **,M(lj ,*Tilt",T 

ent "Enter 2nd Pixel |*,u, "2nd Pixel Scan line #*,M(2) 

cll 'pix' (T/2,P,A,N,R) jell 'GP'(A, •dist'(N,D),Lt3J,L(lJ) 

cll 'pixdir '(A,L(1] ,B) ;L (4] -fasn(sin (B)sin (D) /cos (L( 3] ) ) *L(2] 

cll 'equa'(L(l),L(2] ,S,0,M,L) 

cll 'latlon'(S4-'second'(M(l]-M[2] ) ,L,L[1] ,L(2]) 
cll 'pix'(T/2,Q,A,N,R) 

cll 'pixlat'CLdJ , 'dist'(N,D), •pixdir'(A,L(ll ,B) ,L[5J) 

cll ‘pixlon'(L(5],L(21,D,B,L(6]) 

wtb P,10,10;wrt P,*PIRST PIXEL LAT* , 'dm'( L [3] ) 

wrt P, "FIRST PIXEL LONG", 'da* (L(4] ), "FIRST PIXEL #",P 

wrt P, "FIRST PIXEL SCANLINE t * ,M (1] , " SECONL PIXEL *",Q 

wrt F, "SECOND PIXEL SCAN LINE #" ,M (2) , "TILT" ,Tjgsb " * 

wrt P, "SECOND PIXEL LAT* , 'da' ( L ( 5J ), "SECOND PIXEL LONG" , 'da' ( L (61 ) 
gto "START" 

" 9 *. 

"Enter Reference Pixel Lat-Lon, X-Y corrds, station pix Lat-Long,*: 
"Get station pixel X-Y values": 

ent "Enter Ret Pixel lat.*,L(S] ,*Ret Pixel lon.*,L(6] ,*Ret Pixel #*,Q 

ent "Ret Pixel Scan line •*,M(1] , "Enter Tilt",T 

ent "Enter Station Latitude", L( 3) , "Station Longitude* ,L( 4] 

cll 'pix* (T/2,U,A,N,R) jell *GP' (A, 'di«t'(N,D) ,L(5) ,L(11) 

cll ‘pixdir '(A, L(1J ,B) ;L(6) -fasnCsin (B)sin (0)/cos (L (5] ) ) «L(2| 

cll 'equa'(L(ll ,L(2] ,S,0,M,L);rt«M(21:gsb "Iter* 

M (11 (21 -M) *M;prno (M,0) •*M;ptnd (P,0) *P 

wtb P, 10,10: wrt P,*RbPERENCE PIXEL LATITUDE" , 'oa* ( L (51 ) 

wet P,"KEP PIXEL LONGITUDE", *aa'(L(bl) , "REF PIXEL X VALUE" ,u 

wet F,"Rt.F PlAbL I VALUE", M (1 1 , "STATION LATITUDE" , *da‘ ( L ( 31 ) 

wet F,"biATION longitude", ' da'(L(4)) , "TILT", TjgsD " " 

wet F, "STATION PIXEL X VALUE" , P, "STAT ION PIXEL Y VALUE" ,M 
yto "UTAhT" 


" 10 ": 

"Entee u.P. Lat-Lon, Tilt, Pixel*; eet Pixel Lat-Lon": 


200: ent "Enter G.P. Lat" #L[1J ,"G.P. Lon" ,L(2) , "Tilt" "Pixel# #P 
201: cll •pix*(T/2,P,A,N);cll 'disf(N,D) 


202 

203 

204 
20 5 
206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 
22 3 
224 
22 5 
226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 
24 2 

243 

244 

245 

246 

247 

248 


cll 'pixdir* {A,L(1J ,B);cll 'pixlat • ( L [1] L131 ) 
cll •pixlon*(L(3J »L(2] ,D,B, L[4] ) 

wrt F^,"g!p? latitude" , •dm* ( L(11 ), "G.P. LONGITUDE" , ’dm* (L 12) )» "TILT" ,T 

wrt F, "PIXEL i",P;cll ' • 

wrt F, "NADIR ANGLE" ,N," AZIMUTH ANGLE" , A 

wrt F, "DISTANCE TO PIXEL" , D , "BEARING TO PIXEL ,B 

wrt F, "PIXEL LATITUDE", •dm'(L(3)), "PIXEL LONGITUDE" , ’amM L (4 ] ) 

gto "START" 

" 11 ": 

"Change Roll, Pitcn, and/or Yaw ; 

ent "Enter Roll" ,"Pitcn",X , "Yaw" ,Y 

fmt 2/,25x,"NEW ROLL, PITCH, YAW:";wrt F 

tmt C35," - ",t7.2;wrt F, " ROLL" ,"PITCH" ,X , "YAW" , Y 

gsb " " 

gto "START" 

"Notes" : 

if F«701jfrot 3/;wrt 701 
if F*6;wtb 6,102,125 

«t THE C2NAV PROGRAM CAN BE USED TO CALCULATE A VARIETY OF" 

F,'“CZC£ NAVIGATION PROBLEMS. THE ROUTINES ARE ACCESSED BY USE OF" 
wrt F, "THE SPECIAL FUNCTION KEYS. NOTE: ALL ANGLES ARE ENTERED 

n n»y decimal degrees, output is IN DEGREES MINUTES." 

FUNCTION KEY FO WILL CONVERT THE" 


wrt F, "FOR DEGREES -MINUTES ENTRY, 

wrt f! "ANGLE AND ENTER IT; USE KEY FO INSTEAD OF THE CONTINUE KEY. 
wet F, "TILT IS INITIALIZED TO 20; ROLL, PITCH, YAW TO 0. TILT IS 
wrt F, "ASKED EACH TIME, USE Fll FOR R.P.Y. VARIABLES ARE GLOBAL, 
wtb F, 10,10 

fmt "KEY",30x,"DESCRIPnON",/;wrt F 

fmt 0,"F",f2.0," : ",c;fmt l,7x,c;fxd 1 . 

wrt F, 6, "FROM EsiOATORIAL CROSSING LONG, TIME, AND STATION LAT-LONG, 

wrt F+.l, "CALCULATE TIME OF STATION PASSAGE AND PIXEL 

wrt F, 7, "FROM PIXEL LAT-LONG AND PIXEL t> CALCULATE G.P. LAT-LONG. 

wrt F, B,"FROM 1ST PIXEL LAT-LONG, X-Y VALUES, AND 2ND PIXEL X-Y VALUES, 

wrt F+.l, "CALCULATE 2ND PIXEL LAT-LONG." 

wrt F, 9, "FROM REF PIXEL LAT-LONG, X-Y VALUES, AND STA PIXEL LAT-LONG, 

wrt F+.l, "CALCULATE STATION PIXEL X-Y VALUES." ^ 

wrt F, 10, "FROM G.P. LAT-LONG AND PIXEL I, CALCULATE PIXEL LAT-LONG. 

wrt F, 11, "ENTER NEW ROLL, PITCH, AND YAW VALUES." 

wrt F,5,"DEG MIN - DEC DEG AND TIME - SCAN LINE CONVERSIONS. 

wrt F,0, "DEGREES MINUTES ENTRY." 

if F*701;wtb 701,12 

ret 


15 


249: 

250: 

251: 

252: 

253: 

254: 

255: 

256: 

257: 

258: 

259: 

260: 

261: 

262: 

263: 

264: 

265: 

266: 

267: 

268: 

269: 

270: 

271: 

272: 

273: 

274: 

275: 

276: 


”Var iables” : 

*A » satellite azimuth to pixel”: 

”B « true bearing from G.P. to pixel”: 

”C ■ ”: 

”D ■ distance to pixel in degrees great circle arc”: 

”E ■ ” : 

”F * output device number”: 

”G ■ ” : 

”H ■ ” : 

”I ■ general counter”: 

■J ■ ” : 

”K ■ equatorial longitude used in 'Iter'”: 

”L > equatorial longitude of satellite crossing”: 

”M ■ scan lines from equator to G.P.”: 

”N ■ satellite nadir angle from G.P. to pixel”: 

”0 > degrees great circle arc from equa. to G.P.": 

”P * pixel number (primary or reference pixel)”: 

”Q « pixel number (secondary pixel) ” : 

”R ■ satellite mirror shaft angle": 

”S * seconds since equatorial crossing”: 

”T > satellite tilt angle (2 tau)”: 

”U ■ used in 'Iter'”: 

”V ■ used in ' Iter ' ” : 

”W “ satellite roll angle”: 

”X ■ satellite pitch angle”: 

”y ■ satellite yaw angle”: 

”Z * Time of equatorial crossing”: 

”L(*1 * iati tude- longitude array (odd«lat, even*lon) ” 
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CZSUB.OOC - Documentation of the CZCS Navigation Subroutines 
CZSUBI - CZCS NAVIGATION SUBR #1 

Given a CZCS G.P. latitude, longitude and travel time in seconds, 
this routine calculates a new G.P. latitude and longitude. 

CALL CZSUBI (XINCLP, PERIOD ,XLATO ,XLONGO ,T ,XLAT ,XLONG, I ERR) 


CZSUB2 - CZCS Navigation Subt #2 
Given tne CZCS G.P. latitude, longitude, tilt, and scan angle, 
tnis subroutine will return the pixel latitude and longituoe. 

CALL CZSUB2(XINCLP,ALTC0R,XLAT0,XL0NG0,TILT,XNU,XLAT, 

* XLONG, PITCH, ROLL, XAW, I ERR) 


************************************************************** 

CZSUBI - CZCS Navigation Subr #3 

This subroutine returns tne Latitude ana Longitude or the Grouno 
point of the Satellite. The input variables are the Latitude ana 
Longitude, Tilt and Scan Angle ot a Pixel, ana tne Roll, Pitcn, 
and Xaw, Satellite altitude factor, and equatorial tracKline 
inclination angle ot tne Satellite. See CZSUB Variable inaex for 
explanation of variable names. 

CALL CZSUB3( XINCLP, ALTCOR,XLAT ,XLONG ,TILT,XNU,XLATO , 

* XLONGO, PITCH, ROLL, XAW, I ERR) 


CZSUB4 - CZCS Navigation Subr »4 
Given a satellite G.P. latituae and longituoe, 

this subroutine calculates tne longituae ot equatorial crossing, 
and tne satellite travel time from tne equator to the G.P. 

CALL CZSUB4 (XINCLP, PERIOD ,XLAfU ,XLONGO ,XLONGP ,T, lERR) 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


##**#**^******************************************^*********^******* 

CZSUB5 - CZCS Navigation Subr #5 
Given a pixel latitude and longitude and an equatorial crossing 
longitude, this subroutine iterates to find a satellite G.P. 
latitude and longitude which (1) contains the pixel in its scan 
and (2) crossed the equator at the given longitude. Also returned 
are the scan angle (XNU) and transit time (T) . If the 
fails, meaning no ground point meets the criteria, a value of —9999 
is returned for XNU. 


CALL CZSUB5 (XINCLP, ALTCOR, PERIOD ,XLONGP ,XLAT ,XLONG,TILT, 

* XLATO ,XLONGO ,T ,XNU, PITCH ,ROLL, YAW ,I ERR) 

♦ ♦♦*<i******* ********************************************* *********** 


#***ik******* ****** ************ 




CZSUBA - CZCS Navigation subr A 

This is an internal subroutine called by several of the other CZCS 
navigation subroutines. Given the tilt, scan angle, roll, pitch, 
and yaw, this routine calculates the trackline relatlive bearing 
to the pixel, the effective scan angle, and the earth surface 
distance (in radians) from the satellite G.P. to the pixel. 


SUBROUTINE CZSU.tJA (ALTCOR , TILT ,XNU,RPSI ,RTHS ,RTH , PITCH , ROLL, YAW 
^ ,IERR) 






CZSUB - VARIABLES 

These are the variables named in the call lists in the various CZCS 

navigation subroutines. All variables are REAL*8 except lERR. 

ALICOR The altituae correction (R+h)/R, where R is the radius of 
tne earth and h is the orbital altitude of the satellite. 
Supplied by tne main or calling program. 

lERR This is a logical variable which when true enables error 

trace printouts at tne terminal. Supplied by the main or 
calling program. LOGICAL* 2 

PERIOD The orbital period of the satellite in seconds. Supplieo 
by tne main or calling program. 

PITCH The state of rotation of the satellite about tne Y axis 
in degrees. Supplied by the main or calling program. 

ROLL The state of rotation of the satellite about the X axis 

in degrees. Supplieo by the main cr calling program. 
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RPSI The trackline bearing to the pixel in radians* Returned 

by CZSUBA. 

RTH The earth surface distance from the G.P* to the pixel in 

radians* Returned by CZSUBA. 

RTHS The effective scan angle (as modified by tilt) in radians* 

Returned by CZSUBA* 

T The orbital travel time of the satellite between two 

points in seconds* For example the time for satellite 
transit from the equator to the current G.P* Supplied or 
returned depending on subroutine. 

TILT The tilt angle of the CZCS scan view in degrees. This 

is as read from the satellite (or double the tilt angle 
of the mirror shaft ,Tau) . Supplied by the main or 
calling program* 

XINCLP The trackline equatorial inclination angle, 9.28 deg 
west of north* 

XLAT The latitude of the pixel in decimal degrees* Supplied or 

returned* Or, the new latitude in C2SUB1; returned* 

XLATO The latitude of the G.P. in decimal aegrees. Supplied or 
returned* Or, the input latitude in CZSUBl* 

XLONG The longitude of the pixel in decimal aegrees. Supplied 
or returned* The return longitude in CZSUBl* 

XLONGO The longitude of the G.P* in decimal degrees. Supplied 
or returned* The input longitude in CZSUBl* 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 


XLONGP The longitude of equatorial crossing in decimal aegrees. 
supplied or returned* 

XNU The mirror scan angle in degrees. Supplied by main cr 

calling program* 

YAM The state of rotation of the satellite about the Z axis 

in degrees* Supplied by the main cr calling program. 

Note that while the angular variables are supplied or returned 
in degrees, the internal computer computations are performed 
in radians* These angular variables are converted in the 
subroutines and then reconverted before return. The exception 
to this is CZSUBA whose call list contains radians, variables 
which start with an X are in degrees and those with an R 
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c 

c 


are in radians. For example within CZSUB3 , XNU is converted to 
RNU and XL‘^T becomes RLAT. 


*** LIST OP SUBPROGRAMS 1*1 FILE "C2SUBS" : (C2SUB1 - C2S0BA) 

• CZSUBl 

• CZSUB2 

- CZSUB3 

" CZSUB4 

■ CZSUB5 

" CZSUBA 

*** END DIRECTORY *** 






CZSUBl - CZCS NAVIGATION SUBR #1 

Given a CZCS G.P. latitude, longitude and travel time in seconds, 
this routine calculates a new G.P. latitude and longitude. 


22 July 1981, W. Wilson, visibility Laboratory, 
Scripp's Institution of Oceanography. 


*********************** 






C 


C 


c 


c 

c 

c 

c 

c 

c 


c 


SUBROUTINE CZSUBl (XI NCLP , PERIOD ,XLATO ,XLONGO ,T ,XLAT ,XLONG,IERR) 


IMPLICIT REAL*8(A-H,0-Z) 
REAL*8 PI ,PID2 ,DEGRD,DATAN , 
REAL*8 RINCLP,RLATO,RINCLO, 
LOGICAL I ERR 


DATAN2 ,DSIN ,DCOS ,DARCOS,DARSIN, DSQRT 
DT,X ,RLAM,DTAN, DABS, RLAT 


PI « 4.0*DATAN (1. ODO) 

DEGRD « PI/180.0 

RINCLP * XINCLP*DEGRD /* equatorial inclination 

RLATO - XLATO*DEGRD 

Calculate tracXline inclination at initial G.P. 

RINCLO « DASIN(DSIN(RINCLP)/DCOS(RLATO) ) 

Calculate degrees trackline traveled in T (fixec earth) 
DT ■ T*360. 0/PERIOD*DEGRD 

Calculate new G.P. latitude 

X • DC0S(DT)*DS1N(RLAT0)+DSIN(DT)*DC0S(RLA70)*DC0S(RINCL0) 

RLAT - DASIN(X) 

XLAT » RLAT/DEGRD 
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C Calculate delta longitude 

RIAM - OASIN(DSIN(RlNCLO)/OCOS(IU.Ar)*OSIN(DT)) 

OLAH ■ RLAM/OEGRl) 

New longitude « old long del Ion -f earth rotation 
XLONG « XLONGO OLAM -t- T*4.167E-3 

IP(IERR) PRINT 2000, XLATO XLONGO ,XLAT , XLONG, OLAM, T 
FORMATC CZSUBl - XLATO , XLONGO ,XLAT, XLONG, OLAM, T > ',/, 

* 1P6G13.6) 

RETURN 

*** END SUBROUTINE CZSUBl *** 

END 


CZSUB2 - CZCS Navigation Subr #2 
Given the CZCS G.P. latitude, longitude, tilt, and scan angle, 
this subroutine will return the pixel latitude and longitude. 

22 July 1981, W. Wilson, visibility Laboratory, 

Scripp's Institution of Oceanography. 


SUBROUTIN E CZSUB2 ( XI NCLP,ALlCOR,XLATO ,XLONGO ,TILT ,XNU ,XLAT, 

* XLONG, PITCH, ROLL, YAW, lERR) 

IMPLICIT REAL*8(A-H,0-Z) 

R£AL*8 PI ,Pi 02,DEGRD 

REAL* 8 RINCLP,RLATO ,RINCLO ,RPSIP,RPSI ,RTH ,X, RLAT , OLAM 
REAL* 6 RTHS 

LOGICAL I ERR 

PI « 4. 0* OAT AN (1.000) 

DEGRD « PI/180.0 

IF(IERR) PRINT 2000 
2000 FORMATC ENTER CZSUB2' ) 

Get Psi and Theta from CZSUBA 

CALL CZSUBA (ALTCOR,TILT,XNU , RPSI ,RTHS , RTH , PI TCri ,ROLL , YAW,I ERR) 

RINCLP • XlNCLP*OEGRO 
RLATO ■ XLATO *UEGR0 

Calculate trackline inclination at G.P. latituae 
RINCLO ■ DASIn(DSIN(RINCLP)/DCOS (RLATO)) 

C Calculate true bearing to pixel 

RPSIP ■ RPSI - RINCLO 
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C Calculate pixel latitude 

X ■ OSIH (RLATO)*DCOS(RTH) + DCOS(RIATO) *DSIN (RTH) *DCOS (RPSIP) 
KLKt • DASIM(X) 

XLXr » RIAT/DEGRD /* convert back to degrees 

C 

C Calculate delta longitude 

X ■ DSIH(RPSIP)/DCOS(RLAT)*DSIN(RTH) 

DLM - DASIM(X) 

XOLAM > DLAM/OEGRD 

C Calculate pixel longitude 

XLONG > XLOMGO > XOLAM 

C 

PSIP ■ RPSIP/DEGRD /* convert back to degrees 

TH « RTU/OEGRO 

IP(IERR) PRINT 2003, PSIP, XLAT,XDLAM 
2003 FORMAT (• PSIP ,XLAT ,XDLAM ■ •,1P4G13.6) 

C 

RETURN 

C *** END SUBROUTINE C2SUB2 *** 

END 

..******^.**************.****** 

SUBROUTINE CZSUB3 

This subroutine returns the Latitude and Longitude of the Ground 
Point of the Satellite. The input variables are the Latiti^e and 
Longitude, Tilt and Scan Angle of a Pixel, and the Roll, Pitch, 
and Taw, Satellite altitude factor, and equatorial tracklinc 
inclination angle of the Satellite. 

22 July 1981, W. wilson. Visibility Laboratory, 

Scripp's Institution of Oceanography. 


SUBROUTINE C2SUB3( XINCLP,ALTCOR,XLAT , XLONG ,TILT,XNU,XLATO , 

* XLONGO,PITCH,ROLL,TAW,IERR) 

C 

IMPLICIT REAL*8(A-H,0-2) 

REAL* 8 PI ,PID2,0EGRD 

REAL*8 RPSI ,RTH,RLAT,RINCLP,A,B,C,0,AB,XD,X,XP,X1,X2 
REAL* 8 RLATO 1 ,RLAT02 , RLATC , RI NCLO ,RPS I P ,DLAM 
REAL* 8 RTHS 
LOGICAL lERR 
C 

PI - 4.0*OATAN (l.ODO) 

OEGRO ■ PI/180.0 /* Degrees-Radians conv 

C 

IF(IERR) PRINT 2000 
rOOO FORMAT (* ENTER CZSUB3 ' ) 

CALL CZSUB A ( ALTCOR , T I LT , XN U ,RPS I , RTHS , RTH , P I TCH , ROLL , Y AW , I E RR ) 
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RLAT ■ XLAT*DEGRO 
RINCLP • XINCLP*OEGRD 


Set up Quadratic variables from CZCS angular coeponents, 
See docuaentation for derivation. 

a * Srn** ~ OSIN (RTH) *DS1N (RPSI)*DSIN (RINCLP) 

C - dsin!oth)*dcos(rpsi) 

0 a 000S(RINCLP)**2 


AB a A*B /• r 

XD a b**2 + C**2 

IF(IERR) PRINT 2001 ,A, B, C ,D,XD,AB 
X a AB**2 - XD*{A**2-D*C**2) 

XP a 0.0 

IP(X.GT.l.E-lO) XP a DSQRT(X) 
PORMATC A,B,C,D,XD,AB a MP6G13.6) 


/* Calculate intermediate vals 


X2 Z lAB * vBwiS values from the 

X2 a (AB - XP)/xo /* Quadratic. 

IF(IERR) PRINT 2005,X,X1,X2 
FORMATC X,X1^2-MP3G13.6) 

RLATOl a OASIN (XI) 

RLAT02 a OASIN (X2) 

IF(IERR) PRINT 2006 , RLATOl ,RLAT02 
FORMATC RLATOl ,RLAT02-MP2G13. 6) 

tilt is positive choose the negative quad, 
component and if negative the positive. 

RLATO a RLAT02 

IF(TILT.LT.O.O) RLATO a RLATOl 
XLATOl a RLATO 1/DEGRD 
XLAT02 a PLAT02/DEGRD 
XLATO a RLATO/Dr.GRD 

Calculate true bearing to pixel 
RINCLO a OASIN (OSIN (RINCLP)/DCOS (RLATO) ) 

RPSIP a RPSI - RINCLO 
IF(IERR) print 2002,RPSIP,RTH,RLAT 
FORMATC RPSIP,RTH,RLAT-MP3G13.6) 

Calculate longitude difference G.P. to Pix 
X a OSIN (RPSIP) *0SIN(RTH)/DC0S(RLAT) 

OLAM a OASIN (X) 

XLAM a OLAM/OEGRO 

XLONGO a XLONG ♦ XLAM /* calc G.P. Longitude 

»XLAT01,XLAT02, XLONGO, XLAM 
FORMATC LAT1,LAT2,L0NG,DL0NG - ',1P4G13.6) 
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FSIP - RPSIP/OEGIO) 

IF(IERP) PRINT 2004,PSIP,XLATO,XLAM 
2004 FORMATC PSIP,XLAT,XOLAM - MP4G13.6) 
RETURN 

C *** END SUBROUTINE CZSUB3 *** 

END 


CZSUB4 - CZCS Navigation Subr *4 
Given a satellite G.P. latitude and longitude, 
this subroutine calculates the longitude of equatorial crossing, 
and the satellite travel tine from the equator to the G.P. 

22 July 1981, W. Wilson, visibility Laboratory, 

Scripp's Institution of Oceanography. 


SUBROUTINE CZSUB4 (XINCLP, PERIOD ,XLATO ,XLONGO ,XLONGP ,T, lERR) 
IMPLICIT REAL*8(A-H,0-Z) 

REAL*8 PI ,PI02,0EGR0,0ATAN ,0ATAN2,DSIN ,OCOS, CARGOS, OARSIN,OSURT 
REAL* 8 RI NCLP , RLATO ,RLAM , RT , OTAN , DABS , O'! 

LOGICAL I ERR 

PI • 4.0*DATAN(1.0D0) 

DEGRD « PI/180.0 

RINCLP - X1MCLP*0EGRD 
RLATO > XLATO*DEGRD 

Calculate the delta longitude 
RLAM > DASIN (OTAN (RLATO) *DTAM (RINCLP) ) 

Calculate the degrees traveled and transit tine 
RT - DASIN (OSIN (RLATO )/OCOS (RINCLP) ) 
or > RT /DEGRD 
T ■ DT* PERIOD/ 360.0 

Calculate the equatorial crossing longituoe 
OLAM > RLAM/OEGRO 

XLONGP • XLONGO - OLAM - T*4.l67E-3 
C 

If(IERR) PRINT 2000, XLATO, XLONGO, XLONGP ,0'i , OLAM, T 
2000 FORMAT(' CZSUB4 - XLATO,XLONGO ,XLCNGP , D1 ,DLAM ,T - *,/, 

* 1P6G13.&) 

C 

RE TO RN 
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c 

c$ 


*** END SUBROUTINE CZSUB4 


END 


CZSUBS - CZCS Navigation Subc #5 
Given a pixel latitude and longitude and an 
longitude, this subroutine iterates to find a .*^**^* 

latitude and longitude which (1) contains the pixel . 

and (2) crossed the equator at the given longitude. Also retu^ed 
are the scan angle (XNU) and transit time (T) . It 
fails, meaning no ground point meets the criteria, a value of 
INTL(-9998) is returned for XNU, 

22 July 1981 R w. wilson. Visibility Laboratory, 

Scripp's Institution of Oceanography, 


C 

C 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


L 


SUBROUTINE CZSUBS ( XI NCLP,ALTCOR, PERIOD, XLONGP ,XLAT ,XLCNG,TILT , 
• XLATO,XLONGO,T,XNU, PITCH, ROLL, YAW, I ERR) 

IMPLICIT REAL*8(A-H,0-2) 

LOGICAL lERR 

, Q^O /* initial pixel limits 

XH ■ 1969. 

ITER ■ 0 

Start with a known pixel lat-lon 


10 CONTINUE 

XM ■ (XL+XH)/2.0 

XNU ■ -39.36 ♦ (XM-1.0)*0.04 


/• estimate pixel number 
/• and scan angle 


Go calculate a coresponding G.P. lat-lon 


CALL CZSUB3 ( Xl NCLP, ALTCOR,XLAT ,XL0NG ,TI L i ,XNU , XLATO , 

* XLOnGO, PITCH, ROLL, YAW, 1 ERR) 

Go tinu tne estimateo G.P. 's equatorial crossing Ion 


CALL C^SUb4(XI NCLP , PL f. 100 , X LA TO ,XLUNG0 ,XL0 NGg ,T ,I ERK) 
Compare estimateo eq. Ion with input eq. Ion 


UlP - aLONGO - XLONGP 

IF(DABS(U1P) .Li’. .0002 b) GO TO 100 /» good - go return 

HER • ITEF >1 . __ 

IF ( I'lLR.Gl . 20) GO TO 200 /* wcn t wotR - quit 

or 


25 


on o n nnnnnnnnnnnno 


c 


IF(OIF.GT.O.O) XH * XM 
IF(OIF.LT.O.O) XL • XM 


/* xdju«t acxn limits 
/* for nsxt tty 


2000 


C 

c 


IFtlERR) PRINT 2000 ,XM,XNU ,XLATO ,XLONGO ,XLONGP «XLONGQ, OIF 
FCRKATC CZSUBS • XM,XNU,XLATOKXLONGO ^LONGP ,XLONGQ,OIF 
/,F7.1,F8.2»1P6G11.S) 

GO TO 10 /• go try again 


100 CONTINUE 

IF(IERR) PRINT 2000,XM,XNU,XLATO ,XLONGO,XLONGP,XLONGQ,OIF 
RETURN 

200 CONTINUE 

PRINT 2001 

2001 FORMAT (' FAILED TO CONVERGED - PROBABLY OUTSIDE IMAGE') 

XNU ■ -9999. 

RETURN 

C **• END SUBROUTINE CZSUBS *** 

END 

C$ 


CZSUBA - CZCS Navigation subr A 

This is an internal subroutine called by several of the other CZCS 
navigation subroutines. Given the tilt, scan angle, toll, pitch, 
and yaw, this routine calculates the trackline relatlive bearing 
to the pixel, the effective scan angle, and the earth surface 
distance (in radians) from the satellite G.P. to the pixel. 

22 July 1981, W. Milson and J. Nolten, Visibility Laboratory, 
Scripp's Institution ot Oceanography. 


SUBROUTINE CZSUBA (AI TCOR,TILT,XNU, RPSI ,RTHS ,RTH , PITCH , ROLL, YAW 
• ,IERR) 

IMPLICIT REAL*8(A-H,0-Z) 

REAL*8 PI ,PI02,OEGRD,IROOT 

REAL” 8 PI3D2,PPSI ,RTH,RNU,RTHS,RTILT 

LOGICAL I ERR 


PID2 • DATAN(1.000)*2.0 

PI ■ PID2*2.0 /* calculate pi 

PI3D2 - P1D2*3.0 

OECRO • PI/180.0 /* degree - radian conversion 

IF(IERR) PRINT 20C0 ,ALTCOR , TILT , XNU 
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2000 rOBMATC ENTER CZ5UBA , ALTC0R,T1LT,NU • *«3r8.3) 

RNU ■ XNU*0EGRD 
T • TlLT/2.0 
RT • T*OEGRO 
IP(TILT.NE.O.O) GO TO 1 
RPSI ■ PI02 

IP(»U.LT.O.O) RPSI • PI302 
RTHS - DABS (RNU) 

GO TO 2 
1 CONTINUE 

IF(IERR) PRINT 30Qi,RT,RNU 
30U1 FORMAT (' 1',' RT,RNU MP3G13.6) 

Calculate Thcta**n and Phi-n 

IROOT - 1.000/ OSQRT (2.000) 

RTHR ■ 2.*OACOS(IROOT*(DSIN (RT)*OCOS(RNU)+DCOS(RT))) 

RPHIM - OATAN(OSIN(RNU)/(DSIN(RT)-OCOS(RT)*OCOS(RNU)))+PI 
IF(IERR) PRINT 3002 ,RTHM,RPHIM,RTHR 

3002 FORMATC 2',' RTHM,RPHIM,RTHR MP3G13.6) 

Calculate Theta and Phi 

RTUP • 1. 0O0/DSQRT( 2.000) 

RPHIP - OATAN (DTAN (RTHR)*OSIN (RPHIM) ) 

IF(TILT.LT.O. 0) RPHIP - PI + RPHIP 
IP(IERR) PRINT 3003, RTHP, RPHIP 

3003 FORMAT(* 3',' RTHP.RPHIP* ,1P2G13.6) 

IF (ROLL.NE.O. .OR. PITCH. NE.O. .OR. YAh.NE.O.) GO TO 3 
RPSI a RPHIP 
RTHS a RTHP 
GO TO 2 
3 CONTINUE 


Correct for coll pitch and yaw 

U a OSIN (RTHP )*DCOS (RPHIP) 

V a OSIN (RTHP) *OSIN (RPHIP) 
a a OCOS(RTHP) 
y a YAW'OEGRD 
R a ROLL*OEGRD 
P a P1TCH*DEGRD 

xp a (1.000 ♦ y*R*P)*u ♦ (-y+R*p)*v p*w 
yp a Y*o ♦ V - R*W 

2P a (-P+yaR)*o ♦ (-pay+R)av + ^ 

IF(IERR) PRINT 3004 , U , V, W ,XP ,y P ,ZP 
3004 FORMAT(' 4 ',' U,V,W ,XP,yP,ZP* , 1P6G13. 6) 
RTHS a OACOS(ZP) 

RPSI a OATAN2(yP,XP) 
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2 


2001 

C 

C 

CS 


CONTINUE 

RTH • OASIN (ALTCOR*OSIN (RTHS)) - RTHS 
THS - RTHS/DEGRO 
TH • RTH/OEGRD 
PSI - RPSI / OEGRO 
IF(IERR) PRINT 2001, THS ,TH, PSI 
PORMATC THS, TH, PSI ■ MPlGll.S) 

RETURN 

*** END SUBROUTINE CZSUBA *** 

END 


***CZCS NAVIGATION PROGRAM*** 

PROGRAM CZNAV 

This program uses the CZCS Navigation Subroutines to perform 
* variety of CZCS location calculations. A menu of functions 
is provided, variables are global so that, if iterating, only 
changing parameters need be reentered. Output is displayed 
on the terminal and written to an output file 'CZCS.OATA*. 

22 July 1981, J. Nolten, Visibility Laboratory, Scripp's 
Institution of Oceanography 

IMPLICIT R£AL*8 (A-G,N-Z) 

INTEGER*4 PTR, TTY , CHOICE 
CHARACTER* 80 CHAR 
LOGICAL lERR 

DATA lERR/. FALSE. /,PTR/ 16/, TTY/ 1/ 

/* SETUP 

PER ■ 104.09 /* SATELLITE CONSTANTS 

PERIOD • PER*60. 

XINCLP • 9.28 
C ALTCOR • (952. + 6371.)/ 6371. 

ALTCOR - 1.1494 
REV - 1./. 12375 
ROLL ■ 0. 

PITCH • 0. 

YAW - 0. 

TILT • 20. 

ZERO - 0. 

NADIR - 0. 32 
C 

OPEN (UNIT - PTR, 

* FILE - 'CZCS. DATA', 

* STATUS ■ 'UNKNOWN', 


2t 


c 

1000 

1010 


1020 

1030 

1040 


2000 

2010 

2020 


C 

C 

C 

10 

1050 

1055 

1060 


1070 


ACCESS « 'SBQUENTIAl.' , 

FORM • 'UNFORMATTED', 

ERR » 99) 

FORMAT (///,15X, 'N«lcoa« to the CZCS Navigation pcograa.',/) 

FORMAT (5X,'This pcograa provides a variety of CZCS calculations', 

* /, 'selectable frcai the nenu to follow. Besides the terninal,', 

* /, 'input and output are written to file "CZCS. DATA". ' ) 

FORMAT (5X,'Note. ROLL, PITCH, and lAw are initialized to zero, ', 

* /,13X,'and TILT is initialized to +20.',/) 

FORMAT (5X, 'Please enter Date and/or any information you woulo',/, 
' 'lil^e included as header on your output file. (80 char, avail) ') 
FORMAT (A) 

WRITE (TTY, 1000) 

WRITE (TTY, 1010) 

WRITE (TTY, 1020) 

WRITE (TTY, 1030) 

READ (TTY, 1040) CHAR 

FORMAT (//,15X, 'CZCS Navigation Calculations',/) 

FORMAT (2X,A,//) 

FORMAT (5X,50('-') ,/) 

WRITE (PTR,2000) 

WRITE (PTR,2010) CHAR 
WRITE (PTR,1020) 

WRITE (PTR,2020) 

/* MAIN LOOP 


CONTINUE 

FORMAT (//,20X,'MENU' ,/) 

FORMAT (' Angle Conventions: entry and output in OOO.MMmm' ,/ , 

' North and West positive',/) 

FORMAT (' 1; Enter Pixel Lat-Lon , x-Y coords, get G.P. ', 

'Lat-Lon, Equator data.',/, 

' 2: Enter G.P, Lat-Lon, Pixel #; get Pixel Lat-Lon.',/, 

3: Enter Pixel Lat-Lon, X-Y coords, second Pixel X-Y,', 

/,8X,'get second Pixel Lat-Lon.',/, 

• 4: Enter Equatorial crossing Longitude, Time, ', 

'Station Lat-Lon ;',/,8X, 

'get Station Pixel values and Station Time.',/, 

5* Enter Station Lat-Lon, Calculate Equtorial crossinq*. 

• window. ',/, 

' 6: Enter Pixel Lat-Lon, x-Y coords. Station Pixel ', 

Lat-Lon ,' ,/,8X, 'get Station Pixel X-Y coords.',/, 

DegMin-DecDeg and Vime-Scan Lines conversions', 

/,' 8: Change Roll, Pitch, and Yaw settings.', 

/,’ 9: QUIT'.) 

FORMAT (/, 'PLEASE ENTER NUMBER OF CALCULATION TO PERFORM.') 

WRITE (TTY, 1050) 

WRITE (TTY, 1055) 
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WRITE (rrY,I060) 

WRITE (TTY, 1070) 

READ (TTY,*) CHOICE 
C 

00 TO (20,30,40,50,60,70,80,90,100) , CHOICE 
WRITE (TTY,*) 'INVALID CALCULATION TYPE* 

GO TO 10 
C 

4000 FORMAT (//,' INPUT') 

4001 FORMAT ( 20X, 10 ('-'),/, ' OUTPUT') 

4002 FORMAT (/ ,10X,30 ('-'),/ ) . , 

4010 FORMAT (' Pixel Latitude, Longitude: ' ,2F10. 4) 

4011 FORMAT (' First Pixel: ,5X , 'Latitude and Longitude:', 

* 2F10.4,/,5X, 'X and Y values :' ,216) 

4013 FORMAT (' Tilt «',F7.1) 

4020 FORMAT (' Pixel Number : ' , I 5 , ' , Tilt:',F5.1) 

4030 FORMAT (' G.P. Latitude, Longitude :' ,2F10.4) 

4040 FORMAT (' Equatorial crossing Longitude: ' ,F10. 4,/, 

* ' Travel time since equator: ,F7.1,' seconds.') ^ 

4050 FORMAT (' Equatorial crossing Longitude: ' ,P10. 4) 

4051 FORMAT (' Equatorial crossing time:',F7.2) 

4052 FORMAT (' Station not seen by this Satellite pass.') 

4053 FORMAT (' Station Latitude and Longitude: ' ,2F10. 4) 

4054 FORMAT (' Optimal equatorial window from',F9.4,' to', 

4055 FORMAT*(' Station Pixel X and Y values :' ,216) 

4060 FORMAT (' Second Pixel Latitude, Longitude: ' ,2F10. 4) 

4061 FORMAT (' Travel time between scans :', F7. 1 ,' seconds.') 

4062 FORMAT (' Second Pixel X and Y val ues : ' ,21 6) 

4070 FORMAT (' Time of Station passage: ' ,F8. 3,' , Pixel #',I5) 

4090 FORMAT (' Roll, Pitcn, Yaw changed to:',3F8.3) 

5000 FORMAT () 

C . 

C From Pixel into, calculate G.P. into 

C 

20 CONTINUE , 

300U FORMAT ( 'Enter Pixel Latitude, Longitude, Pixel Number, and Tilt ) 

WRITE (TTY ,3000) 

READ (TTY,*) XLAT ,XLONG, PI XNOl ,TILT 
XLAT > OMTOD(XLAT) 

XLONG > OMTOO(XLONG) 

XNU - -39.04 * .04*PIXNOl 

CALL C2SUB3( XI NCLP,AL1C0R, XLAT, XLONG, TILT, XNU, XLATO, 

* XLONGO, PITCH, ROLL, YAW, I ERR) 

CALL C2SUB4 (XINCLP , PERIOD , XLATO , XLONGO ,XLONGP ,T, lERR) 

XLAT • Ol'ODM(XLAT) 

XLONG - DTODM (XLONG) 

XLATO ■ DTODM (XLATO) 

XLONGO - DTODM (XLCNGO) 

XLONGP > DTODM (XLONGP) 
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I 

I 

? 

! 


WRITE (TTY, 4 000) 

WRITE (TTY ,4010) XLAT,XLONG 
WRITE (TTY ,4020) PIXN01,TILT 
WRITE (TTY, 4001) 

WRITE (TTY,4030) XLATO,XLONGO 
WRITE (TTY,4040) XLONGP ,T 
WRITE (TTY ,4 002) 

C 

WRITE (PTR,4000) 

WRITE (PTR,4010) XLAT,XLONG 
WRITE (PTR,4020) PIXN01,TILT 
WRITE (PTR,4001) 

WRITE (PTR,4030) XLATO ,XLONGO 
WRIl-E (PTR,4040) XLONGP ,T 
WRITE (PTR,4002) 

C 

WRITE (TTY,*) 'Use RETURN Icey to continue pcogran when ready.' 
READ (TTY, 5000) 

GO TO 10 
C 

C From G.P. Lat*-Lon and Pix I, Calculate Pixel Lat-Lon 

C 

30 CONTINUE 

3010 FORMAT ('Enter G.P. Latitude, Longitude, Pixel «, and Tilt') 
WRITE (TTY, 3010) 

READ (TTY,*) XLATO ,XLONGO ,PIXN01 , TILT 
XLATO * OMTOO (XLATO) 

XLONGO > OMTOO (XLONGO) 

XNU « -39.04 * .04*PJ.XNOl 

CALL CZSUB2(XINCLP,ALTCOR, XLATO, XLONGO , TILT, XNU, XLAT, 

* XLONG, PITCH, ROLL, YAK, I ERR) 

XLAT « OT-OOM(XLAT) 

XLONG • OTOOH (XLONG) 

XLONGO • OTOOM( XLONGO) 

XLATO • OT-OOM (XLATO) 

WRITE (TTY, 4 000) 

WRIT-E (TTY, 4030) XLATO , XLONGO 
WRITE (TTY, 4020) PIXN01,TILT 
WRITE (TTY, 4001) 

WRITE (TTY, 4 010) XLAT, XLONG 
WRITE (TTY, 400 2) 

C 

WRITE (PTR,4030) XLATO , XLONGO 
WRITE (PTR,4020) P1XN01,TILT 
WRITE (PTR,4001) 

WRITE (PTR,4010) XLAT, XLONG 
WRITE (PTR,4002) 

C 

WRITE (TTY,*) 'Use RETURN Key to continue program wnen ready.' 



31 


oooo 


READ (TTY, 5000) 
GO TO 10 


Prom Lat-Lon and X-Y for Pixel #1 and X-Y tor Pixel *2 
Calculate Lat-Lon for Pixel #2 

40 CONTINUE . , ^ 

3020 FORMAT ('Enter Pixel #1 Latitude# Longitude# Pixel ai^ Scan #s') 

3021 FORMAT ('Enter Pixel and Scan Is for Pixel 12, and Tilt') 

MRITE (TTY, 3020) 

READ (TTY,*) XLAT,XLONG, PIXNOl ,SCAN01 
XLAT - OMTOO(XLAT) 

XLONG « OMTOO (XLONG) 

WRITE (TTY, 3021) 

READ (TTY,*) PIXN02 ,SCAN02 ,TILT 
XNO » -39.04 + .04*PIXNOl 
T » (SCAN01-SCAN02) *0.12375 

CALL C2SUB3(XINCLP,ALTCOR,XLAT,XLONG,TILT,XNU,XLATO, 

* XLONGO,PIi-CH,ROLL,YAW,I£RR) 

WRITE (TTY, 4000) 

WRITE (TTY,4011)DTODM(XLAT) ,DTODM( XLONG) , P I XNO 1 ,SCANO 1 
WRITE (PTR,4000) 

WRITE (PTR,4011) DTODM(XLAT) ,DTODM(XLONG) ,PI XN01,SCAN01 
CALL C2S0B1 (XINCLP , PERIOD ,XLATO,XLONGO ,T , XLAT , XLONG, lERR) 


C 


C 


C 


XLATO > XLAT 

XLONGO - XLONG 

mi « -39.04 + .04*PIXNO2 

CALL CZSUB2(XINCLP, ALTCOR, XLATO, XLONGO, TILT, XNU,XLAT, 
* XLONG, PITCH, ROLL, YAW, I ERR) 


XLAT « OTODM(XLAT) 
XLONG > OTOOM (XLONG) 
XLATO - DTODM( XLATO) 
XLONGO > O'l'OOM (XLONGO) 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 


(TTY, 4062) 
(TTY, 4013) 
(rrY,4001) 
(TTY, 4060) 
(TTY, 4061) 
(TTY, 4002) 
(PTR,4062) 
(PTR,4013) 
(PTR,4001) 
(PTR,4060) 
(PTR,4061) 
(PTR,4002) 


PIXN02,SCAN02 

TILT 

XLAT, XLONG 
T 

PIXN02,SCAN02 

TILT 

XLAT, XLONG 
T 


WRITE (TTY,*) 'Use RETURN Key to continue program when ready.' 
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READ (TTY, 5000) 
GO TO 10 


Given Eq crossing Ion and Sta. Lat-Lon, 
Calc time of sta passage, scan posit 


SO CONTINUE 

3050 fORKAT ('Enter Equatorial Crossing Longitude, Time (HH.MNb)*) 

3051 FORMAT ('Enter Station Latitude, Longitude, and Tilt') 

WRITE (TTY, 3050) 

READ (TTY,*) XLCNGP.EQT 
XLONGP « OMTOO(XLONGP) 

WRITE (TTY, 3051) 

READ (Tl'Y,*) XLAT , ALONG, TILT 
XLAT • OMTOD(XLAT) 

ALONG « OMTOD (ALONG) 

CALL C2SUB5(XINCLP,ALTCOR, PERIOD, XLONGP, XLAT, ALONG, TILT, 

* XLATO ,XLONGO ,T ,XNU , PITCH ,ROLL, YAW ,I ERR) 


C 


C 


STT‘ • DTODM(DMTOD (EOT) + DMIOD (T/6 000. ) ) 
XLONGP « DTODM( XLONGP) 

XLAT « OK)ON(XLAT) 

ALONG > OTODM (ALONG) 

PIANOS « (XNU * 39. 04)/. 04 

WRITE (TTY, 4000) 

WRITE (TTY, 4053) XLAT, ALONG 
WRITE (TTY, 4050) XLONGP 
WRITE (TTY, 4051) EQT 
WRITE (TTY, 4 001) 

IP (INTL(XNO) .EU.-9998.) THEN 
WRITE (TTY, 4052) 

GO TO 51 
END IF 

WRITE (TTY, 4070) STT, PIANOS 

51 CONTINUE 
WRITE (TTY, 4 00 2) 

WRITE (PTR,4000) 

WRITE (PTR,4053) XLAT, ALONG 
WRITE (PTR,4050) XLONGP 
WRITE (PTR,40S1) EQT 
WRITE (PTR,4001) 

IF (INTL(XNU) .EQ.-9998. ) THEN 
WRITE (PTR,4052) 

GO TO 52 
END IF 

WRITE (PTR,4070) STT, PIANOS 

52 CONTINUE 
WRITE (PTR,4002) 
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WRITE (TTY,*) 'Use RETURN key to continue program when ready.' 
READ (TTY, 5000) 

GO TO 10 


Prom Station Latitude, Longitude, Calc Equatorial window 


C 


c 


c 


c 

c 

c 


60 CONTINUE 

WRITE (TTY, 3051) 

READ (TTY,*) XLATO ,XLONGO ,TILT 
XLATO > DMTOO( XLATO) 

X1.0NG0 > OMTOD(XLONOO) 

CALL C2SUB4 (XINCLP, PERIOD ,XLATO ,XLONGO ,XLONGP ,T, lERR) 
CALL C2SUB2( XINCLP, ALTCOR, ZERO ,XL0NGP, TILT, 35. 0D0,XLAT, 

* XLONGE , PITCH, ROLL, YAW,IERR) 

CALL CZSUB2(XINCLP ,ALTCOR, ZERO ,XL0NGP, TILT, “35. 0D0,XLAT, 

* XLONGW,PITCH,ROLL, YAW,IERR) 

XLONGP a DTODM(XLONGP) 

XLONGE > OTODM (XLONGE) 

XLONGW - OTODM (XLONGW) 

XLAT « DTODM(XLAT) 

XLATO « OTODM (XLATO) 

XLONGO - OTODM (XLONGO) 


WRITE (TTY ,4000) 
WRITE (TTY ,4053) 
WRITE (TTY ,4001) 
WRITE (TTY ,4040) 
WRITE (TTY ,4054) 
WRITE (TTY ,4002) 


XLATO ,XLONGO 

XLONGP ,T 
XLONGE ,XLONGW 


WRITE (PTR,4000) 
WRITE (PTR.4053) 
WRITE (PTR,4001) 
WRITE (PTR,4040) 
WRITE (PTR,4054) 
WRITE (PTR,4002) 


XLATO ,XLONGO 

XLONGP ,T 
XLONGE ,XLCNGW 


WRITE (TTY,*) 'Use RETURN key to continue program when ready.' 
READ (TTY, 5000) 

GO TO 10 


From reference pixel and station Lat-Lon get sta. X-Y 


70 CONTINUE 

WRITE (TTY ,3020) 

READ (TTY,*) XLAT ,XLONG, PIXNOl ,SCANO 1 
XLAT > l)MTOD(XLAT) 

XLONG - DMIOD(XLONG) 

WRITE (TTY, 3051) 

READ (TTY,*) SLAT ,SLONG,T ILT 
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SLAT > CMTOD(SLAT) 

SLONG > 3MT00(SL0NG) 

XNU • -39.04 ♦ ,04*PIXN01 

CALL CZSUB3 ( XI NCLP ,ALTCOR,XLAT ,XLONG ,TI LT ,XN0 , XLATO , 

* XLONGO,PITCU,HOLL, YAM.IERR) 

CALL C2SOB4 (XINCLP, PERIOD ^LATO ,XLONQO ^LONGP ,T, lERR) 

T1 ■ T 

CALL CZSUB5(XINCLP,ALTCOR,PERIOD,XLONGP, SLAT ,SLONG, TILT, 

* XLATO ,XLONGO ,T ,XNU, PITCH ,ROLL, YAW ,IERR) 

PIXK02 - (XNU + 39. 04)/. 04 

SCAN02 ■ SCANOl - (T - Tl)/0. 12375 
C 

XL AT ■ DTODM(XLAT) 

XLONG <• OTODM(XLONG) 

SLAT - DTODM(SLAT) 

SLONG - DTOOM (SLONG) 

C 

WRITE (TTY, 4000) 

WRITE (TTY, 4011) XLAT,XLONG,PIXN01, SCANOl 
WRITE (TTY,4C53) SLAT, SLONG 
WRITE (TTY,<i013) TILT 
WRITE (TTY, 4001) 

IF (INTL(XNO) .EQ.-9998) THEN 
WRITE (TTY, 4052) 

ELSE 

WRITE (TTY, 4055) PIXN02,SCAN02 
END IP 

WRITE (TTY, 4002) 

C 

)«RITE (PTR,4000) 

WRITE (PTR,4011) XLAT , XLONG, P I XNOl, SCANOl 
mRITE (PTR,4053) SLAT, SLONG 
WRITE (PTR,4013) TILT 
WRITE (PTR,4001) 

IF (INTL(XNU) .EQ.-9998) THEN 
WRITE (PrR,40S2) 

ELSE 

WRITE (PTR,4055) P I XN02 ,SCAN02 
END IF 

WRITE (PTR,4002) 

C 

C 

WRITE (TTY,*) 'Use RETURN key to continue pcogram when ready.' 
READ (TrY,5000) 

GO TO 10 
C 

C DegMin-DecDeg and Time-Scan Lines conversions 

C 

80 CONTINUE 

3080 FORMAT ('Enter value to be converted, space, and identifier:') 
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o n o o 


Seconds ' , 

Deciaal Oeg cees ' , 


3083 

3084 

3081 

3082 

81 


3085 FORMAT (T15,'* L" ■ Scan lines* ,/ ,T15 , S" 

* /,T15,'" DM" ■ Degrees Minutes' ,/ ,T15 , DD 

* /,T15,’ "Q- » QUIT*) 

FORMAT (F7.1»* seconds equals' »I5f' scan lines. ) 

FORMAT (I5»' scan lines equals' ,F7.1 , ' seconds.') . , , 

FORMAT (F10.4,' degrees .minutes e ,oals' ,P10. 4, 

FORMAT (F10.4,' decimal degrees equals' ,F10. 4, ' degrees. minutes ) 
IFIRST «‘l 
CONTINUE 
WRITE (TTY, 3080 

IF (IFIRST. EU.l) WRITE (TTY, 3085) 

READ (TTY, 1040) CHAR 
IF (INDEX(CHAR, 'Q') .NE.O) GO TO 10 
READ (CHAR ,* ) VAL 
IF (INDEX(CHAR, 'M') .NE.O) THEN 

WRITE (TTY, 3081) VAL ,DMT0D(VAL) 
t SE IF (INDEX (CHAR, 'S' ) .NE.O) THEN 
,.riTE (TTY, 3083) VAL, VAL/0. 12375 
ELSE IF iINDEX(CHAR, * L' ) .NE.O) THEN 
WRITE (TTY, 3084) V AL, VAL* 0. 12 375 
ELSE 

WRITE (TTY, 3082) VAL, DTODM ( VAL) 

END IF 

WRITE (TTY, 4 00 2) 

IFIRST ■ IfIRST + 1 
GO TO 81 


Change Roll, Pitch, Yaw values 


VU UOWlinwc . «r.oTu»r npropcciM 

3090 FORMAT ('Enter new Roll, Pitch, Yaw vaiues. DECIMAL DEGREES) ) 
WRITE (TTY, 3090) 

READ (TTY,*) ROLL, PITCH, YAW 
WRITE (TTY, 4090) ROLL, PITCH ,YAw 
WRITE (PTR,4090) ROLL, PITCH , YAW 
WRITE (PTR,4002) 

C 

GO TO 10 
C 

99 CONTINUE 

WRITE (TTY,*) 'ERROR ON FILE OPENING 
STOP 'EXIT C2.CS Navigation' 

C 

100 CONTINUE 

ENOFILE (PTR) 

CLOSE (PTR) 

STOP 'EXIT CZCS Navigation' 

FND 
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REAL FUNCTION OMTOO*8(X) 

REAL*8 X,P3,P2,OINT,OHOO,OSIGN,OABS 
REAL*b ONTOO 

P3 ■ OABS(X) 

P2 « OINT(P3) + OMOO(P3,1.0O0) *100.0/60.0 
U4TOD « P2*OSlGN (1.0D0,X) 

RETURN 

END 

REAL FUNCTION OTOOM*8(X) 

IMPLICIT REAL*8(A-U,0>Z) 

REAL*8 P3,P2,OINT,ONOO,X,OSIGN,OABS 
P3« OABS(X) •»- .0000005 

P2 ■ DINT(P3) ♦ DMOD(P3,i.OOO) *60. 0/100.0 
OTOOM « P2*OSIGN (1.000, X) 

RETURN 

END 


